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Abstract 

In this paper we analyse the vacuum fluctuations of quantum fermionic current densities, 
(j A (x)}, induced by a magnetic flux running along the core of the conical space produced 
by the presence of cosmic string, in a (1 + d)— dimensional bulk, for 2 < d < 5. In order 
to develop this investigation we construct an expression for the fermionic propagator, which 
can be applied for those different dimensions of the spacetime. Two specific analysis for 
the induced current densities are presented: (i) for massless fields in three and four dimen- 
sions, and (u) for massive fields in all dimensions considered. In the latter specific values 
for the parameters which codify the presence of the cosmic string and the fractional part of 
the ratio of the magnetic flux by the quantum one are considered. Although being a very 
special situation and somewhat unnatural, the corresponding analysis of induced fermionic 
current densities under this circumstance may shed light on the qualitative behavior of these 
quantities in a more general situation. 

PACS numbers: 98.80.C*g, ll.10.G7i, 11.27. + d 



1 Introduction 



It is well known that different types of topological defects may have been formed in the early 
Universe after Planck time by a vacuum phase transition [H [2] . Depending on the topology of 
the vacuum manifold these are domain walls, strings, monopoles and textures. Among them, 
cosmic strings and monopoles seem to be the best candidates to be observed. Cosmic strings are 
linear topologically stable gravitational defects. The gravitational field produced by a cosmic 
string may be approximated by a planar angle deficit in the two-dimensional sub-space. The 
simplest theoretical model describing an idealized cosmic string, i.e., straight and infinitely 
thin, is given by a delta-type distribution for the energy-momentum tensor along the linear 
defect. This object can also be described by classical field theory. By coupling the energy- 
momentum tensor associated with the Maxwell-Higgs system investigated by Nielsen and Olesen 
in [3] with the Einstein equations, Garfinkle [4] found static cylindrically symmetric solutions 
representing vortices, as in flat space-time, and shown that asymptotically the two dimensional 
spaces perpendicular to the vortices correspond to a plan minus a wedge. Moreover, their core 
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have a non-zero thickness with a non-vanishing magnetic fields inside. In additional publication, 
Linet [5] obtained, as a limiting case, exact solutions for fields equations, corresponding to the 
spacetime produced by an idealized cosmic string, with the conical parameter being related with 
the energy per unity length of the vortex. 

Though the recent observation data on the cosmic microwaves background have ruled out 
cosmic strings as the primary source for primordial density perturbation, they are still candidate 
for the generation of a number of interesting physical effects such as gamma ray burst [6], 
gravitational waves [7] and high energy cosmic ray [8]. Moreover, in the framework of brane 
inflation [9j-[llJ, cosmic strings have attracted renewed interest, partly because a variant of their 
formation mechanism is proposed. 

Although topological defects have been first considered in four-dimensional spacetimes, they 
have been analysed in higher-dimensional spacetimes in the context of braneworld. By this 
scenario the defects live in a n— dimensions submanifold embedded in a (4 + n)— dimensional 
Universe. The cosmic string case, corresponding to a two additional dimensions, has been 
analysed [121 13j. There it was shown that the gravitational effects of global strings can be 
responsible for compactification from six to four spacetime dimensions, naturally producing the 
observed hierarchy between electroweak and gravitational forces. 

The analysis of quantum vacuum fluctuations associated with matter fields in curved space- 
time are relevant in a semiclassical theory of gravity p3|. The validity conditions under which 
the vacuum expectation value (VEV) of the energy-momentum tensor can act as the source for a 
semiclassical gravitational field are discussed in |15U16j . In this context, the VEV of the energy- 
momentum tensor in the idealized cosmic string spacetime, have been considered in |17j-|21j and 
|22|-[25j. for scalar and fermionic fields, respectively. It has been shown that these quantities 
depend on the parameter which codify the conical structure of the two- geometry. Moreover, 
considering the presence of magnetic flux running along the cosmic strings, there appears an ad- 
ditional contributions to the corresponding vacuum polarization effects associated with charged 
fields [26]-[3T]0 Also by imposing boundary conditions on quantum fields, additional shifts in 
the vacuum expectation values of physical quantities, such as energy density and stress take 
place. This is the well-known Casimir effects. The analysis of Casimir effects in the idealized 
cosmic string space-time have been developed for a scalar [33 1. vector [33] and fermionic fields 
|35j . obeying boundary conditions on the cylindrical surfaceso 

Additionally to the contribution for the energy-momentum tensor associated with charged 
quantum fields, the presence of a magnetic flux along the cosmic string induces vacuum current 
densities, as well. This phenomenon has been investigated for massless scalar field in [28] and 
recently for massive one in [29] . There the authors have shown that a vacuum induced azimuthal 
current densities, (j^), take place if the ratio of the magnetic flux by the quantum one has a 
non- vanishing fractional part. H So in a more general analysis about the influence of vacuum 
polarization effects on the semiclassical theory of gravity, induced currents produce additional 
contribution. This quantity is among the most important characteristics of the vacuum state. 
Though the corresponding operator is local, due to the global nature of the vacuum, the VEV of 
the current density describes the global properties of the bulk and carry important informations 
about the structure of the defect core and the magnetic flux. In addition to describe the 
physical structure of the quantum fields at a given point, the current density acts as source of 
the electromagnetic field in the Maxwell's equations. It therefore plays an important role in 
modeling a self-consistent dynamics involving the electromagnetic fields. 

Recently the fluxes by gauge fields play an important role in higher dimensional models including braneworld 
scenario (see for example [32]). 

2 Also vacuum polarization effects induced by a composite topological defect has been analysed in |36| . 

3 The analysis of electric current densities induced by magnetic flux have been investigated in [37]- [ID] in a flat 
spacetime. 
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Although the analysis of the VEVs of the energy-momentum tensor associated with charged 
fermionic fields in a cosmic string spacetime having a magnetic flux line running along its core, 
has been developed by many authors; as far as we know the investigation of the corresponding 
induced current densities has been restricted for a three-dimensional spacetime only in |41j . So 
the objective of the paper is to extend this analysis as general as possible. Adopting this line 
of investigation, here we shall calculate the vacuum fluctuations of fermionic current densities 
in a (1 + d)— dimensional cosmic string spacetime, considering massless and massive fields. As 
we shall see these current densities depend crucially on the fractional part of the ratio of the 
magnetic flux by the quantum one. 

This paper is organized as follows: In Section[2]we briefly present some important results ob- 
tained in [42j, associated with fermionic propagators in an arbitrary (1 + d)— dimensional cosmic 
string space-time, with d > 2, in the presence of a magnetic flux. Although the main objective 
is to investigate induced currents in a higher-dimensional spaces, the formalism obtained allows 
us to develop similar analysis for a three-dimensional spacetime, as well. Moreover, in this 
section we also present a closed expression for massive fermionic Green function, expressed in 
terms of finite sum of Macdonald function, K Vl for specific values of the parameters which codify 
the presence of cosmic string and the fractional part of the ratio of the magnetic flux by the 
quantum one. In Section [3] we calculate induced vacuum fermionic current densities for mass- 
less and massive fields. For massless fields, we explicitly calculate all components of fermionic 
current densities for three and four dimensions and shown that only azimuthal components are 
different from zero. For massive fields and for specific conditions above mentioned, we calculate 
also all components of fermionic currents for d = 2, 3, 4, 5. We shown that only for d = 2 
there appears a non-vanishing charge density. In Section 0J we summarize the most important 
results obtained. The Appendix [A] contains some details related with the analysis of radial and 
azimuthal current densities associated with massless fields in three-dimensional spacetime. In 
this paper we shall use h = G = c = 1, denote the electric charge by e and use for the metric 
tensor the signature +2. 

2 The fermionic propagator 

In this section we briefly review the obtainment of the fermionic propagator derived in previous 
paper [42] . in the (1 + d)— dimensional cosmic string spacetime in the presence of a magnetic flux 
running along its core. By using the generalized cylindrical coordinate with the cosmic string 
on the subspace defined by r = 0, being r > the polar coordinate, the corresponding metric 
tensor is defined by the line element below: 

ds 2 = g MN dx M dx N = -dt 2 + dr 2 + a 2 r 2 dtp 2 + ^{dx l f . (1) 

% 

The coordinates system reads: x M = (t,r,ip,x l ), with cp G [0, 2ir], and t, x l G (— oo, oo). The 
parameter a, smaller than unity, codify the presence of the string and is given, in a four dimen- 
sional spacetime, in terms of the linear mass density of the string by a = 1 — 4/i. Specifically for 
a six-dimensional space, i = 3, the geometry represented by ([1]) corresponds to two dimensional 
conical surface transverse to a 3-flat brane. 

In order to calculate the fermionic propagator, we shall start the construction of the for- 
malism admitting that the spacetime is six-dimensional. So we shall adopt for the 8x8 Dirac 
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matrices, the representation below, constructed in terms of the 4x4 ones 

1°\ yr= ( f • 7 \ rv = _}_( ■ 7 
7° j ' \r-i j ' ar\0-i 

7 D' r -C?). r -(-o)' < 2 » 

where 75 = Z7°7 1 7 2 7 3 , I represents the 4x4 identity matrix and f and (p stand the ordinary unit 
vectors in cylindrical coordinates. This set of matrices satisfies the Clifford algebra {T M , F N } = 

~ 2 9 I (8)- 

Let us now introduce an infinitesimally thin magnetic flux running along the string. This 
will be implemented by considering a six-dimensional vector potential below 

A M = Ad M <f , (3) 

being A = Jp. 

The fermionic propagator on curved spacetime and in the presence of a gauge fields, can be 
expressed in terms of a bispinor, Vp(x,x'), by [H] 

S F (x,x' ) = (if + e4 + M)V F (x,x') , (4) 

where the covariant derivative operator reads 

J7 = T M (d M + U M ), (5) 

being Um the the spin connection given in terms of the T— matrices 

n M = -\r N v M r N , (6) 

and 

4 = r M A M . (7) 

Moreover, the bispinor must obey the differential equation below [31], 
[□ - ieg MN (D M A N ) + ieY, MN F MN - 2ieg MN A M V 



N 

1, 

I' 



-e 2 g MN A M A N -M 2 -^TZ V F (x,x') = — -^5 6 (x - x')I {8) , (8) 
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with 



E MiV = 1 [T M^ N] Dm = Vm _ ieAM (9) 
1Z represents the scalar curvature and the generalized d'Alembertian is expressed by 

□ = 9 MN V M V N = g MN (d M V N + n M Vjv - {mn}Vs) ■ (10) 

For the system under consideration the differential operator fC, which acts on the left hand 
side of (JSJ, reduces to 

2/6 Adm - S^A 2 - M 2 , (11) 



a 2 r 2 ^ a 2 r 2 
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where 



J (8) 



E 3 
S 3 



with 



a d 
cr 3 



an 



r a z r A v y 



a = -d\ + a 2 + -d r + 



(12) 



(13) 



Because £ 3 g ^ is a diagonal matrix the bispinor, T>p(x,x') can be given in terms of a 2 x 2 

bispinor, T>p'(x, x ), in a diagonal form. The latter obeying the 2x2 matrix differential operator 
below: 



A + 



(1-aV 3 ^ 



-4-^t(1 - af + -^(1 - a)Acr 3 
4a 2 r 2 a 2 r 2 



2ie 



Aft, 



a 2 r 2 



t4 — 



vP(x,x') 



-g 



5\x-x')I {2) 



(14) 



On basis of these results, for this six-dimensional cosmic string space-time, the fermionic 
propagator can be expressed as shown below: 



Sf(x, x') 



iT°d t + iT r d r + iT^d^ + tPa 



1 — a q e<& 
2 ( § ) 2vr 



D F (x,x') . (15) 



The calculation of the induced current can be given in terms of the Euclidean bispinor. 
The latter been related with the ordinary Feynman propagator by the relation T>e(t, r; r', r') = 
—H)f(x,x'), with t = it [H]. In the following we shall consider the Euclidean bispinor only. 

In order to calculate the 2x2 Euclidean bispinor, T> E (x, x'),we shall first obtain the complete 
set of normalized eigenfunctions of the Euclidean version, IC, of the differential operator given 
in ([E 



£$ A (x) = -A 2 $ A (x) 



with A 2 > 0. In this way V^\x, x') is given by 
V^ { x,x') = Y J ^ {X) ^ [X ' ) 



X 2 



A 2 



ds^2$\(x)<f>\(x') e 
x 2 



(16) 



(17) 



The eigenfunctions of (1161) can be specified by a complete set of quantum number associated 
with operators that commute with JC and among themselves. They are: p T = —id T , pi = —idi 
for i = 3, 4, 5 , L v = — id^ and the spin operator, 173. Let us denote these quantum numbers 
by (k T , k l , n, a), where (k T , k l ) € (—00, 00), n = 0, ±1, ±2, ... , a = ±1, and p, that 
satisfies the relation A 2 = p 2 + k 2 + M 2 , assuming value in the interval [0, 00). Imposing on the 
eigenfunctions the Dirichlet boundary condition on the stringH 



$ A (i,r = 0,cp,x) = , 
the corresponding eigenfunctions are: 

Ak.x 



(18) 



$1 } (a0 



[a(2^) JV + 1 ] 1 /2 



Ak.x 



[a(2vr) 7V + 1 ] 1 A 



(+) 



a (pr)uj 



(-) 



(19) 



4 For this geometry the only non- vanishing spin connection is II ^ = J(l — a)S 3 



(8)- 

A more general boundary condition implies singular solutions in r = 0. It is our intention to extend the 
analisys presented in this paper admitting singular solutions. 
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where 

^ +, = (J)' wH = (?) ' (20) 

In the above expressions represents the Bessel function, and = n ± ^— 2 ~~ 7' where we 
have defined e^4 = ej^ = 7, being 7 £ (0, 1). 

Finally with the help of [15] we can express T> E (x,x') in a integral representation: 

Vl/«(™7 2 *) \ 

V |/a (rr'/2 S )) ^ ' ^ 

where I^{z) is the modified Bessel function. 

Although we have developed this formalism for a six-dimensional cosmic string space-time, 
we can adapt it for a three, four and five dimensional spaces. The reason resides in the repre- 
sentations adopted for the Dirac matrices. For three dimensional space the corresponding Dirac 
matrices are given in terms of 2 x 2 Pauli matrices; and for a four and five dimensions, they 
are given in terms of the 4x4 off-diagonal matrices of (|2|). Moreover, in (|2ip we should assume 
N = 1, 2, 3, 4, which are related with the dimensions of the space considered. 

Unfortunately, it is not possible to integrate over the variable s in (|2ip and to express the 
Euclidean bispinor in a closed form. In fact all the calculations of vacuum polarizations con- 
sidering massive fields involve very complicated integrals. So in order to present more workable 
expressions, the massless limit of the corresponding results are adopted. Here in this paper, 
we shall explicitly calculate the renormalized induced fermionic current for different dimensions 
of the space-time in closed expressions. So, in order to do that we shall consider two different 
limiting cases which allow us to reach this objective as shown in the following subsections. 

2.1 Massless fermionic fields 

In this subsection we provide the bispinor associated with a charged massless fermionic fields in 
an arbitrary dimension. Taking M = in (|2ip we obtain [45], 



Vf{x,x') 



1 1 



a(2ir) 2" (rr') 2 (— sinhw) 2 



£ e *K*-V)| %l\l*-i/2^ u ) ^° I (22) 



where 







(Ax) 2 + r 2 + r' 2 



Q\v-\/a-l/2( COShu )) 



cosh-u = ■ . (23) 

2rr' y ' 

In the above expression, is the associated Legendre function. 

For a three dimensional space, N = 1, the Legendre function can be expressed in a integral 
form: 

Q u _ 1/2 (caBku) = -= / dt == . (24) 

v2Ju v cosh t — cosh u 

Substituting the above expression into (f22j) . it is possible to obtain a closed expression for the 
bispinor by developing the sum over the quantum number n. Admitting that the parameters 
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S + = — 2 °' > — 7 and 5 = — + 7 that appear in v±, are both restricted to the interval, 
< 5^ < 1, the Euclidean bispinor reads: 

£ V 7 47T 2 a v / 277A yWit-coshu V S^>(t) J v ; 

where 

5 (±) (t) = g^!2 smh(^t/g) - sinh[(^ - l)t/g] _ 
cosh(i/a) — cos(c^ — 92') 

For a four dimensional space, N = 2, the Legendre function can be represented by a simpler 
expression: 

Ql /2 i /2 ( coshu ) = *y| J^r^ ■ ( 2? ) 

In this case, substituting this function into ()22[) and developing similar procedure as before, we 
obtain 

£>IV,zO = 8 ^2 arr / sinhn ( ( } 5(-)(„)) ) ' (28) 
where S , ^(u) has the same functional form as S , ^ ± ^(i) in 
2.2 Massive fermionic fields 

In the analysis of vacuum polarization effects associated with quantum scalar fields in a cosmic 
string spacetime, Smith [18], Davies and Sahni |19j . and Souradeep and Sahni [20], have shown 
that when the parameter a is equal to the inverse of an integer number, i.e., when a = ~, 
being q an integer number, the corresponding Green functions can be expressed in terms of 
q images of the Minkowski spacetime functions. Recently the image method was also used in 
|33] to provide closed expressions for massive scalar Green functions for a higher-dimensional 
cosmic string space-time. The reason for the application of this method is because the order of 
the Bessel functions which appear in the derivations of the Green functions becomes an integer 
number. Unfortunately, for the fermionic case, the order of the Bessel function depends on the 
factor — coming from the spin connection. However, if we consider a charged fermionic field 
in the presence of a magnetic flux running along the string, an additional factor will be present, 
the ratio of the magnetic flux by the quantum one, 7. For the case where 7 is equal to — 
and a = 1/q, the order of the Bessel function becomes an integer number, and in this case 
it is possible to use the image method to obtain the fermionic propagator in a closed formJll 
Although being a very special situation and somewhat unnatural, once that by a typical grand 
unified theory a realistic value for a is very close the unity, the analysis of induced fermionic 
current under this circumstance may shed light on the qualitative behavior of this quantity for 
non-integer q. 

In order to derive a more compact expression for the Euclidean bispinor in the above men- 
tioned situation, let us go back to the integrand of (|2ip . the heat kernel: 

K(x x'- s) - e ~ (A3 ° t;+ " V e m ^ ( VlA*K/2*) \ 



6 The expressions for the bispinor in five and six dimensions are given in [42]; however, because we only analyse 
the three and four dimensional case for massless fields, we shall not repeat them here. 

7 This special situation was also used in [46] to analysis fermionic vacuum polarizations by composite topological 
defect. 
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Defining 

K(+\x,x';s) =J2^ nA(p Iw+\/a(rr'/2s) , (30) 

n 

with u + = n + — 2°' > ~~ 7) the order of the modified Bessel function becomes an integer number. 
For this case, with help of \A7\ H2] we can write: 

(x, x';s) = J2 e inA ^I nq (rr'/2s) = - cos (^ + ^) . (31) 

Analogously, defining 

K(-\x,x';s) =Y,e inAip I\v-\/a(rr'/2s) = e iA ^ e mA ^I nq+1 (rr' /2s) , (32) 

n n 

after some steps we obtain 

K^{x,x';s) = e< 1 - 1 ^^ - ^ co "\~ + ~ ) e~~ . (33) 



k=0 



Finally substituting and ) into the heat kernel and integrating on the variable s we 

find with the help of 



M N / 2 9-1 1 / 1 



l/g)AvP e -2i7rfc/g / i 



fe=0 

with 

p£ = (Ax) 2 + r 2 + r /2 - 2rr cos (Ay>/g + 2vr£:/g) . (35) 

Now having derived the basic expressions to calculate the fermionic propagator, we leave for 
the next section the main objective of our paper, that is the calculations of induced fermionic 
current densities. 



3 Induced current densities 

The fermionic current density operator is an expression bilinear on the field. By using the 
general definition for fermionic propagator, 

iS F {x,x') = (0\T(V(x)V(x'))\0) , (36) 

with * = ^r , we can express the VEV of fermionic current density as the trace of the product 
of the Dirac matrix by the fermionic propagator shown below: 

(j A (x)) = -ie lim Tr{T A S F (x' ,x)) . (37) 

x'— >x 

In the two next subsections, we shall present the steps developed to obtain vacuum expec- 
tation values of the fermionic current densities for massless and massive fields. 



3.1 Massless field 

In this subsection we shall calculate the induced fermionic current densities for a massless fields, 
considering three and four dimensional cosmic string spacetimes. 
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3.1.1 Three dimensional case 



From the general expression for fermionic propagator given by (|15p . for a three-dimensional 
spacetime, it is express in terms of a 2 x 2 bispinor as shown below: 



Sf{x , x) = i 



ia 3 d t - a r d r - —a*d 9 - J-^—^ a <P^ + Jl^ 
ar ar 2 ar 



V%\x',x). (38) 



Let us start the calculation by the time component of the current density. Because the 
bispinor above is in a diagonal form, in the calculation of charge density, only the time 
derivative term survives when we take the trace over the matrices product o^Sf{x' , x). So the 
corresponding VEV is given by 

(j°(x)) = ie lim Tr(d T V^\x', x)) . (39) 

x'— >x 

Here we shall write the bispinor in the form, 
The relevant calculation here is 

d T I^(u)=d T f = S&(t). (41) 

Ju V cosh t — cosn u 

Let us first take Ar = Aip = 0. Introducing a new variable t := 2 arcsinh(y/\/2), the function 
/(^(u) reads: 

I^(z) = 2 I" dV S(±) (2arcsinh(y/V2)) 
7^ + y/y*-{z-l) 

with z = cosh u = 1 + tttt • 

In order to calculate the VEV of the charge density operator, we shall write the above 
function as the sum of two terms: 

!<*>(*) = I?\z) + I?\z) = 2[ 1 dV jEM 



OC 



+ 2 / - y — w , 43 

where S^\y) = S^(t(y)). Subtracting and adding into the integrand of the first term of 
the expansion 

§i±Hy) ~^ + 0(y) (44) 



we get 



where 



V2 + y 2 y 



7 (±) = j&fin + 7 (±) S m S ^ (45) 



A (±)/4 » = 2r * f^4-2) (46) 

7vi=T Vv - ( z - 1) V V 2 + y 2 y / 
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and 



4 



being 77 = arctg ( J 2— 1 1 . In the limit Ar — > 0, the above expression behaves as 



(±)sing 



dy 



a 



2ar\ 



- ^y 2 -(z-l)y 



(47) 



l[ ±)sin9 (z) 



\/2' 



rair 



At 



+ 0(1) 



Substituting this expression into (|40|) . the singular part gives, 



1 



Air At' 



-I 



(2) 



(48) 



(49) 



which coincides with the Euclidean bispinor for a flat three-dimensional space-time. This is 
the only term that contains divergence in the calculation of the charge density. To obtain a 
finite result, we should extract its divergent part. This can be done in a manifest form by 
subtracting from the complete bispinor, the standard Euclidean Green function given above. So 
the renormalized bispinor reads: 



V ( §\x',x) Ren 



1 



/ 1 (+)/in (z)+^'.:. 



r(+), 











I 



(-)/*» 
1 



(50) 



Defining a new variable b = V ' z — 1, we can write: 



(±)/in 



dy 



So, 



6 yjy 2 -b 2 



d T l[ ±)fm (z) 



with /(±)(y) 



W2 + y 2 



1 d 
-V2db 



dy 



f {±) (y) 



Vv 2 



b 2 



(51) 



(52) 



Because the integrand is divergent at the point y = b, we have to change the variable y —> by 
before applying the Leibnitz formula for derivative of an integral. After taken the first order 
derivative, defined x = by and developed some intermediate steps, we find: 



d_ 

db 



dy 



f {±) (y) 



dx d_ ( /W(x) \ _ 

^ x 2 _ b 2 dx [ X J ^/l^l? 



(53) 



Now we can take the limit b — > 0. The term inside the bracket gets a finite result, consequently 
the complete expression vanishes in this limit. The analysis of the time derivative of I^iz) in 
the limit Ar — > is very simple and a vanishing result is promptly obtained. 

As conclusion we can affirm that there is no charge density induced by the magnetic flux for 
the physical situation investigated. 

As to the radial current, < j r >, two distinct contributions are relevant in the product, 
a r Sf(x' ,x). They are expressed as shown below: 



Tr (YS F (x',x)) =Tr 



1 



or 



Vf{x',x) 



(54) 



First we analyse the contribution due to the radial derivative. It is divergent in the coincidence 
limit, and this divergence coincides with the corresponding one given by the standard Green 
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function. Renormalizing this contribution in the same way as we did before we get a finite 
and non- vanishing result H The second contribution is due to the other three terms. The latter 
is finite at the coincidence limit. Taking both contributions together we obtain the imaginary 
expression bellow: 



(fix)) 



le 



Ren. 



dv 



167r 2 r 2 Jo sinh(v/2) cosh(u/a) — 1 



G(v, a, 7) , 



(55) 



with 



4cosh(i>/2a) sinh(t>/2) cosh(try/a:) 2sinh(w) 

H 7(1 — cosh(u/aj) 



cosh(w) 



1 



87 sinh(i>/2a) sinh(i>/2) sinh(t>7/a) 



a- 



(56) 



In appendix [A] it is proved that (|55p vanishes for 7 smaller than | , which is compatible with 
the condition imposed for the parameter 5^ in (|26p . In fact, because we are implicitly assuming 
that a < 1, the condition on the parameter 5 + = (1 — a)/2 — 7, imposes that 1/2 — 7 must 
be bigger than zero. Consequently our analysis is valid only for < 7 < 1/2; moreover, the 
condition 7 < 1/2, naturally appears under specific situation obeyed by the parameters a and 
7 given in the beginning of subsection 12.21 

In the calculation of azimuthal component of the current density we have: 



Ti(^S F (x',x)) 



-Tr 



ar 



a 3 d r + 



i 

ar 



9 v + t;(1 



a )a 



^7 



V%\x',x) 



(57) 



The first contribution due to the radial derivative is given by: 



d r V { *\x',x) 



1 



where 



r(±) 



+ 



(«) 



4n 2 a 
1 



d 



1 



4ir 2 a \J 2rr' dr 



/(+)(«) 
dr \^/2r~?J\ V I^\u) 
1 d ( /(+)(«) 

/(-)(«) 



dt 



v 7 cosh t — cosh u 



S^\t) 



(58) 



(59) 



At this point we may taken Ar = 0, providing coshu = r 2 "tT , and also Atp = 0. 

In the first contribution of (158ft we get a divergent result from the terms inside the matrix 
when we take the coincidence limit. The later provides 



j(±)( 



1 



dt 



-S^(t) 



(60) 



y/2 Jo sinh(t/2) 

In order to obtain a finite result, we should renormalize the bispinor by subtracting from it its 
corresponding Euclidean flat bispinor, 

S(0)/ ' * 1 



f>W(x',x) 



4^Ar /(2) 



(61) 



Expressing the above Green function in a integral form, we get, for the renormalized contribution, 
the following expression: 

^)-<?(0)(t) o \ 

£M_ 5 <o)( t ) J ' 



IF 
16ir 2 r 2 J 



dt 

sinh(t/2) 



(62) 



In order to make a short discussion, some details involved on this calculations are similar to that will appear 
in the next one. 
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where S^(t) is obtained from (|26p by taking 5^ = 0. 

The second contribution of (|58p comes from — gjr~^- Using a procedure similar to the 
adopted in the calculation of the charge density, we can prove that d r I^ L ^ m (u) and d r I^ (u) 
go to zero at the coincidence limit. So, the only effective non-vanishing contribution is given by 

As to the contributions due to the azimuthal derivative, spin connection and magnetic flux in 
(|57p . they provide a finite results in the coincidence limit, so it is not required any renormalization 
procedure. Finally in the coincidence limit we have: 



where 

S (±) (t) = ±5 ± [smh(5 ± t/a) + sinh[(l - <5 ± )t/a]] =p smh(5 ± t/a) . (64) 

Now let us now write down the complete result for the azimuthal current. It reads 

e f°° dt 1 

{f{x))Ren. = Q 2 3 3 / ■ 1 7T7o\ ■ ZTTTTT, s F (^ a ^) » ( 65 ) 
g n 2 a d r d j Q sm h(t/2) smh. z (t/2a) 

where 

F(t, a, 7) = cosh(i/2) [cosh(t/2a) sinh(t7/a) - 27 sinh(t/2a) cosh (£7/0-)] . (66) 

By the above expression we can easily see that F(t, a, 7) vanishes for 7 = 0. The behavior of 
the the azimuthal current as function of the parameters a and 7 can only be provided numerically. 
So in figure[T]we exhibit the behavior of 87r2 J 3 (jf (%)) Ren as function of a and 7. H. 

The appearance of the above induced current is a direct consequence of the well known 
Aharanov-Bohm effect. In principle it should be a periodic function of 7; however, because the 
integrand of (|65p has been obtained under specific conditions, for a fixed value of < a < 1 the 
dependence of < j^{x) >R en with the parameter 7 is positive and convex downwards, presenting 
a maximum value at some specific point. (This fact is also observed by numerical analysis.) 



dt 



8n 2 ar J sinh(t/2) cosh(t/a) 

S(+)(t) 
5(~)(i) 



(63) 



3.1.2 Four dimensional case 

For a four dimensional spacetime, the Dirac matrices in the cosmic string spacetime, can be 
given by 

r = (0) pr = f . - i rV = . - &nd r z = (3) _ ( 6? ) 

ar 

The corresponding fermionic propagator can be obtained by (I15p . being the bispinor the 4x4 
diagonal matrix below 

^ (X ' X)= ( ^ 2) (x',x) J ' (68) 



9 Because the expression found for (|26[) has been obtained admitting that < S ± < 1 our numerical analysis is 
according to this conditions 
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Figure 1: This graph provides the behaviors 87r ^ r3 (jf '(x)) R en , as functions of a and 7, admitting 
these parameters varying in the interval [0.30, 0.50] and [0.05, 0.20], respectively. 

T> E (x',x) being given by 

(x> x) - 1 ( S(+) («) , ? (69) 

with 

(±) gg^V) sin h(^n/«) - sinh[(^ - l)u/a] 
S W = w , s 7 > (70) 

cosn(u/a) — cos(</? — ip ) 
± (1-a) At 2 + Az 2 + r 2 + r /2 

= =F 7 and cosn-u = . (71) 

2 2rr 

Because the bispinor is expressed in terms o elementary functions, the calculations for all 
components of the induced current are relatively simple. Below we present the main steps of 
the calculations. 

Let us first analyse zero component of the current. For the charge density the only non- 
vanishing trace contribution is given in term of time derivative: 

Tr(7 5 F (x', x)) = -Tr(d T V E (x', x)) . (72) 

Here we can take first Acp = Az = Ar = 0, consequently in the coincidence limit, 

d T — — S(±>( U )_ +0 {{t'-t)) . (73) 

smhri [t' — t) 6 

This divergence can be removed by subtracting from the bispinor above the standard Green 
function for a flat four dimensional spacetime. As consequence we obtain a vanishing result for 
the charge density. 

For the radial current we have: 



Tr(-f r S F (x',x)) = iTrj -id r + -j- (^ {3) d v + - a) - i^T^j V E (x',x)^ 



(74) 
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The first contribution is due to the radial derivative. We may take Aip = Az = At = 0, so in 
the coincidence limit we have, 

^* (±) («) — - (737)3 6^T " + Oir - r>) . (75) 

Substituting back this expression into the bispinor, the first term on the right hand side above, 
coincides with the corresponding one given by the standard Green function. Again applying the 
renormalization procedure adopted in the previous calculations, a finite contribution is left. It 
is given below: 

tomdrV E i x ^)Ren. - i8n 2 a 2 r 3 ^ 1 - a 2 + 65" (S- - 1) J ■ [7b) 

As to the second contribution we can analyse the operator below, 

a% + 1(1 - a) - z 7 a 3 ) vf(x>,x) ^ x 
5 + [l + 5 + (28 + -3)] 



5~ [1 + 8- (2d- -3)] J • (77) 

Now considering both contributions above together and taking the trace, we get a vanishing 
result for the radial current. 

As to the azimuthal current, the contributions which provide a non- vanishing traces are given 

by: 

Tv(^S F (x',x)) = -iTr| Y?d r + ±- (^-i(l-a)S 3 + 7 ^) V E (x',x)\ . (78) 

Comparing the above expression with the similar one obtained for the radial current, we can 
see an important difference between them: The presence of the overall matrix — i£ 3 multiplying 
the left hand side of (|74p . Because the steps to be followed are formally the same as the previous 
one, we shall not repeat them. Taking the trace of all contributions, we obtain a real and non- 
vanishing value to the azimuthal current: 

°'" (x)W= l^W (1 + 3a2 - 474) - (79) 

Because the expressions found for (|28|) has been obtained admitting that < 8^ < 1, the above 
result can only be valid for < 7 < 1. 

Finally as to the current along the z— direction, we have: 

Ti( 1 z S F (x' 1 x)) = Ti(d z V E (x',x)) . (80) 

Taking here A</? = Ar = At = 0, we obtain coshu = 1 + tNt, consequently in the coincidence 
limit we may write: 

1 / 1 \ Anr 2 

dz -—S^(u) —> — -3 + 0((z - z')) . (81) 

sum it (z — z') 6 

The first term on the right hand side of the above expression coincides with similar one obtained 
from the standard Green function. So by applying on this calculation the same renormalization 
procedure that we have applied until now, we get ci vanishing result for (j^(^))ijen- 
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3.2 Massive case 

Here we shall analyse the VEV of fermionic current densities associated with massive fields in 
high dimensional cosmic string space times under the circumstance specified in 12.21 In order to 
do that we shall use the bispinor (|34p to calculate the fermionic propagator. However, before to 
start the calculation, it is interesting to analyse this bispinor in the coincidence limit. We can 
verify that taking x' — > x this bispinor is divergent, and that this divergence comes exclusively 
from the k = component. So, in order to obtain a finite and well defined results for the 
currents, we should apply some renormalization prescription. As in the previous calculation, we 
shall apply this procedure in a manifest form, by subtracting from (I34D its k = component. 
So the renormalized bispinor is express as, 

V^(x',x) Ren . = ^ N/2+i ^—^K N/2 (M Pk ) ^ Q ei(1 _ 1/?) A^ e - 2i7 rfc/ g J • (82) 

Adopting TV = 1, 2, ,3, 4 we can construct the fermionic propagator for spacetimes of three, 
four, five and six dimensions. 



3.2.1 Three dimensional case 



The renormalized fermionic propagator in a conical three dimensional spacetime, can be ex- 
pressed by: 



Sf{x ,x) Ren . = i 
being 

V^\x',x) Ren 

with 



ia 3 d t - a r d r a v du 



i (1 — a) 



ar 



ar 



a v a 3 + — ja v + M 
ar 



V^\x',x) Ren . ,(83) 



M i/2 ^ ! ( 1 



e i(l-l/<?)A(/? e -2i7Tfc/<j 



p{ = (At) 2 + r 2 + r' 2 - 2rr' cos (Aip/q + 27rk/q) 



(84) 



(85) 



The first calculation to be realized is for the charge density. We may take, Aip = Ar = in 
the above expressions. The only non-vanishing trace of matrix, ^°Sf(x', x), are given below: 



Tr( 7 °S' i ,(x , ! x)) iie „. = iTr{[id T + Ma 3 ]V^(x',x) Ren } 



(86) 



Let us first to analyse the time derivative of the bispinor. In the coincidence limit, we promptly 
obtain: 



lim d T —K 1/2 (M Pk ) = . 



(87) 



The second contribution is due to the mass term. In the coincidence limit, after some interme- 
diate steps, we find: 



Tr 



iMa 3 vW(x',x) Ren . 



iM 
4?rr 



q-l 



^2 sin(7rA;/g) e 



2 Air sin(irk/q) 



k=l 
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Substituting this expression into (|86p . and using the definition for the VEV of current density 
(f37|) . we obtain: 



eM 



(f(x))Ren = J^^^i^/l) e ~ 



-2Mr sm(nk/q) 



(89) 



k=l 



which presents an exponentially suppressed behavior. 

Now we can calculate the induced electric charge by integration as shown below: 



Q 



drdtp^g^(f(x)) = — (q - 1) = -7 . 



(90) 



The next calculation is the radial current. The non- vanishing contributions to this current 
are given below: 



Tr (YS F (x',x)) =Tr 



d r + — (a 3 d v + U\ -a)- ija :: 



V%\x',x) 



(91) 



For the specific case under consideration, two differential operators naturally emerge. They 



arc: 



D + = 3 r + + id + ) = d r + ^ 



(92) 



because 5 + = 0, and 



D = d r — — (dtp — id 



(93) 



with 5- = 2=1 
q 

The above operators act on the each diagonal components of the bispinor as follows: 



1 



K 1/2 (Mp k 



d r 



1 



K l/2 (M Pk ) 



H o v 

r 



1 



K 1/2 (M Pk ) 



(94) 



and 



K 1/2 (Mp k )e 
-K 1/2 (Mp k )e- 2 ^ 



i(l-l/g)Atp e -2nk/q 



dr 



K 1/2 (M Pk )e 



i(l-l/q)A<p e -2nk/q 



K 1/2 (Mp k )e 



i{l-l/q)A l p e -2 7 rk/q 



(95) 



Because the dependences of the functions inside the brackets, on the radial polar coordinate 
and angular variable, are trough p k , we may use the identity: 



d r = sm(irk/q)- — and d v = - oasiixk/q)- — 
dp k r dp k 

Doing this, we can obtain simplified results for (|94p and (|95|) : 



(96) 



K 1/2 {Mp k ) 



-iwk/q 



d Pk 



K 1/2 (Mp k )e 



i(l-l/q)A<p e -2nk/q 



-ie 



-ink/q 



dpk 



K 1/2 {Mp k ) 



K 1/2 (Mp k ) 



(97) 
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Now returning to the calculation of the radial current density, we have to take the trace of the 
resulting matrix. That means that we have to add both results above, consequently we find 

(j r )Ren = 0. 

For the azimuthal current, the relevant terms to be analysed are: 



Tr (^ Sf (x',x)) = -^Tr i d r + ^ (a% + ~(1 - a) - i 7 a 3 ^ a 3 V^\x',x)\ 



(98) 



In this calculation we can adopt the same approach that we have used in previous calculation; 
however, here there exist an important difference when compared with the later. The matrix <r 3 
appears multiplying the bispinor. For this case, repeating the same steps as we did before, we 
get a non-vanishing result when take the trace of final result. We obtain: 



9-1 



lim Tr OfSr^x)) = V l±^^NA e -*Mr M «k/ q ) . (99) 

x'->x ' 8irr 6 z — ' sm{-irk/q) 

Finally the current density reads, 



,. , sin(7rA;/g) 



= ^ E 1+2M ;fW q h -™r^) . (10 o) 



We may analyse this expression in two limiting cases: 

• In the limit Mr >; 
the leading term is 



In the limit Mr » . , / , , ; , the main contributions come from k = 1 and k = q — 1 and 



{f(x)) Ren « g^e-^sin^) . (1Q1) 

This component presents an exponential suppressed behavior. 
• Massless limit, 

(f{x)) R en = ^3 E -TT7T • ^ 102 ) 

In appendix [X] it is shown that the general expression for the the azimuthal current density 
in three-dimensions for massless fermionic fields, Eq. (|65p . reproduces the result above 
when we substitute in that expression a = 1/q and 7 = (q — l)/2q. 

3.2.2 Four dimensional case 

For a four dimensional spacetime, the renormalized 2x2 bispinor reads, 

V^\x\x) Ren . = ™ g ±- Kl (M Pk ) ( I ei(1 _ 1/9) it e _ 2 ^ /9 ) , (103) 
k=i Hk v y 

with 

p| = Ar 2 + Az 2 + r 2 + r' 2 - 2rr' cos (Ay?/g + 2-irk/q) . (104) 
On the other hand the fermionic propagator is, 

Sf(x, x') = i 



T mi ' 



V E {x,x') , (105) 
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bemgV E (x,x') = diag(V^(x', x), V^(x',x)). 
For the calculation of charge density, we have, 

Tr( 7 5 F (x',x)) &n . = iTv{[id T + M-y°]V(x', x) Ren } . (106) 

Taking the time derivative of the bispinor and then the coincidence limit, we get a zero result. 
As to the last term in above expression we get a vanishing result because the trace of the product 
j°V(x', x)R en , is zero. So as final result we have = 0. 

For the radial current, we can also apply the same procedure as we did in the three di- 
mensional case. The reason is because we can reduce the present analysis for a 2 x 2 matrix 
formalism. Moreover, because we did not have to explicitly calculate the derivative on the 
Macdonald functions in (I97p . we can infer by the same reasons as before, that (j r ) = 0. 

For the azimuthal current density we have: 



Tr(r^5 F (x',x)W = -^Tr| d r + ^ (£ 3 d v + - * 7 X 3 ) T?V E {x,x')\ 



(107) 



Once more, applying to the derivative of the radial and azimuthal coordinates of the bispinor 
the same procedure used before, we finally obtain, after some intermediate steps, the following 
result for the azimuthal current density: 

eqM q ~ ' 



k=l 



MKq (2Mr sm(irk lq) ) - . . Ki (2Mr sin(vr/c lq)) 

r sm(Trk/q) 



(108) 



Two limiting cases deserve to be analyzed: 
• In the limit Mr » . . / , , s , 

sinn(7rK/g) ' 



(TO) - ~^k\ T^T^z- 2Mr Sin(7r/9) • (109) 



2 7r 2 r 2 Y Mr sin(ir /q) 
• Massless limit, 

5-1 

(W) = rh E ■ 2 ( T, \ ■ ( n °) 

The summation on the right hand side of (lllOp can be written in a closed form. Defining 
the sum 

«-i 1 

In(x) =Y J —, , ■ 7 r , (HI) 

sm n (x + t7r/g) 

for even n, I n (x) can be developed by using the formulas |33j 



. , I'J(x) + n 2 I n (x) r . . o 2 1 . 

J n+2 s = nl ; ^ 7 being J 2 (x = -4— - — 2~ — . (112) 
n(n + 1) sm \qx) sin (i) 



Because hi®) = ^ 3 ^ , we have 

{f{x)) = i£h {q2 - l) ■ (113) 

The above result is an analytical function of q, and by analytical continuation is valid for 
all arbitrary values of qf^ 

For the calculation of current along the z— direction we have: 

Tr( 7 3 5 F (x',x)) = Tr(d z V E (x,x')) . (114) 

Taking the derivative of the bispinor and the coincide limit we get a vanishing result. 

10 Substituting a = i and j = 3=1 into JTS} we obtain JTT3J. 
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3.2.3 Higher-dimensional spacetimes 



The main objective of this section is to calculate fermionic current densities in a five and six- 
dimensional cosmic string spacetime. 

In a generalized (1 + d)— dimensional idealized cosmic string spacetime the non vanishing 
singular curvature is located on the (d — 2)— brane, i.e., the scalar curvature has support only 
on the r = 0: R = 47r(a _1 — l)5^ 2 \f). For d > 4 this correspond to a (d — 2)— hyper surf ace. 
Considering a magnetic flux along the brane, which has also support on the brane, an azimuthal 
induced current on the transverse two-dimensional conical surface defined by the coordinates r 
and if according to (pQ), takes place. By using the Maxwell equation, a magnetic field is also 
induced outside the bTane, being directed along it0 Specifically for a six-dimensional bulk, our 
world can be represented by a flat four-dimensional spacetime having a non-vanishing magnetic 
field in it. In this context the analysis which will be developed in this subsection, may be useful 
in a modeling a self-consistent dynamics involving charged quantum fields. 

a an azimuthal electric current will be induced on the transverse two-dimensional cone. 

For a five dimensional case, the Dirac matrices are 4x4 ones and can be represented by, 



r° = 7 (°), r r = r ■ 7 , = — ■ 7 , V x = 7® and T y 



ar 



«75 



(115) 



As to the bispinor, it can be written in a diagonal form in terms of the 2x2 bispinor below, 



^J_if 3/2 (Mp fc ) 
fe=l Pk 



1 

e «(l-l/<3)A</J e -2i7rfc/<j 



with 



Pk 



At 2 + Ax 2 



Ay 



2 i 2 

r 



'2 



2rr' cos (A(p/q + 2-irk/q) 



(116) 



(117) 



The procedure to calculate all components of the fermionic current density here, is very 
similar with the corresponding one for a four-dimensional space-time analyzed in subsection 
13.2.21 To calculate the charge density we have 



Tr(-/°S F (x',x)) Ren . =iTr{[id T + M 1 °}V{x',x) Ren .} . 
Taking first the time derivative of the /c-component of the bispinor, we obtain 



-±-K V2 {Mp k ) 
Pk 



At d 

Pk dp k 



-L K3/2 (M Pk ) 
Pk 



(118) 



(119) 



which goes to zero in the coincidence limit. As to the second contribution of (|118p . the matrix 
j°T>(x', x)_Ren. is traceless. So we conclude that < j°(x) >= 0. 

For the radial component of the current we need to calculate Tr( 7 r S'(x, x')). For this case 
the effective operator which acts on the bispinor is 



1 



8 r + — S^cl 



'1 



a 



ar 



ar 



ar 



(120) 



The contribution due to the extra fifth-coordinate on the Dirac operator is given by ^y r ^d y , 
which provides an off-diagonal matrix. So we conclude that there is no essential difference in 
the calculation of the radial current density when compared with the corresponding calculation 
in four dimension. Here we also obtain < j r (x) >= 0. 



n It our intention to examine the behavior of the induced magnetic field, and its consequence as source of the 
Einstein equations in a backreaction formalism. 
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The components < j x (x) > and < j y (x) > are also both equal to zero. The reasons for 
these facts are similar what happened in the calculation of charge density. The derivative with 
respect to the coordinates x and y on the the bispinor, provide terms proportional to Ax and 
Ay, respectively, which go to zero in the coincidence limit, and the products of 7 X and 7 y by 
T>(x! , x)ji en , are traceless. 

Finally we can see, by using the Dirac operator, and after some intermediate steps, that the 
only non-vanishing component for the fermionic current density is the azimuthal one: 

9^1 -2Mrsin(7rfc/ g ) 

0' V (*)> = ^4-tE -3/ (3 + 6Mrsin(vrfc/g)+4M 2 r 2 sin 2 Wg)) . (121) 

Below we present two limiting cases for the above current: 

• Large mass limit, 

enM 2 p-2Afrsin(7r/ (? ) 

W) « a 2 3 ■ r m • (122) 
47r z r d sm(7r/gj 

Here also there appears an exponential suppressed behavior. 

• Massless limit, 

k=l v ' ' 

For a six dimensional spacetime, the Dirac matrices have already been defined in (JS}. The 
bispinor is also a 8 x 8 diagonal matrix: 

/ V^\x',x) \ 

V^\x',x) 

V%\x',x) 

V Pg } (x',x) / 



T> E (x',x) 



(124) 



being 



with 



p( 2 )(x',x) fen . = ^ Y^^K 2 {M Pk ) ( J e i(l-l/ q )2p e -2i*k/ q ) > ( 125 ) 



pf. = Ar 2 + Ax 2 + Ay 2 + +Az 2 + r 2 + r /2 - 2rr' cos (A<p/g + 27rA;/g) . (126) 



The calculations of all components of the fermionic current in this spacetime will be briefly 
presented in the follows. For the charge density we have, 

Tr (T°S F (x,x')) Ren = iTv{[id T + MT°] V(x',x) Ren .} . (127) 

The time derivative of the /c-component of the bispinor goes to zero in the coincidence limit. 
Moreover, the matrix T°T>(x', x)n en . is traceless. So we conclude that < j (x) >= 0. 

For the radial component we need to calculate Tr (T r Sf(x, x')). The relevant operator is 

- ldr + -IeJ^ + Jl^IzA _ Htf . (128) 
ar y ' ar 2 ar y ' 
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This operator is a 8 x 8 matrix version of the (|120p . Consequently, by considering previous 
analysis, we can conclude that < j r (x) >= 0. 

The components < j l (x) >, for i = 3, 4, 5, given in terms of TV (r i Sp(x, a;')-Ren.) , provide 
vanishing results. The reasons are exactly the same as before: the contributions due to the 
derivative with respect to the coordinates x % on the bispinor go to zero at the coincidence limit 
and the others, given by the product T t T A T>(x', x)n en _, for A / i, and T l T>{x' ', x)R en . are traceless 
matrices. 

Finally we can see that only the azimuthal component of the fermionic current is different 
from zero. After some steps we find: 

M ' ^ 1 

(f(x)) = ^r^V — o- —[2Ki(2Mrsm(irk/q))+2Mrsm(nk/q)K (2Mrsm(7rk/q) 

^ r sua. (irk/ q) 

+ M 2 r 2 sin 2 (vrA;/g))Xi(2Mrsin(7rA;/g))] . (129) 
Finally we present below the two limiting cases: 

• Large mass limit, 

eqM 5 / 2 e - 2Mrsi n(ir/q) 

ir(X)) " 4^/2,7/2 sin 3/2 ( ^ /g) • ( 13 °) 

• Massless limit, 

0'^)> = -rh E ■ 4, ■ / x • (131) 

4tt r D sm (irk/q) 

The above summation can be expressed in terms of an analytical function of q, by ap- 
plying the recurrence relation given in (I112|) taken n = 2 and the expression for l2{x). 
Although In{x) is a long expression, we only need its value for zero argument. It is 

/4(0) = ^ — ^| consequently, 

'" (X))= 185^ (92 - 1)(92 + 11) - (132) 



4 Conclusion and Discussions 

In this paper we have investigated the vacuum expectation value of fermionic current densities 
in a (1 + d)— dimensional cosmic string spacetimes, for 2 < d < 5, induced by the presence 
of a magnetic flux line running along its core. In this analysis two limiting cases have been 
considered. The case where the fermionic fields have no mass, has been considered in three 
and four dimensions only. There we have shown that only the azimuthal components of the 
current are different from zero. Specifically for three dimension, we have founded an integral 
representation for it. In order to provide a better understanding about this current we presented 
its numerical behavior as function of the parameters which codify the presence of the cosmic 
string and the ratio of the magnetic flux by the quantum one. For a four-dimensional spacetime 
we could present the azimuthal component of the current density in a closed form. Considering 
the particular situation where the parameter a is the inverse of an integer number, and the 
ratio of the magnetic flux by the quantum one, 7, is equal to — , massive fermionic fields 
were considered in all dimensions above specified. For this case fermionic propagators could 
be express in terms of finite sum of modified Bessel function, K^, and we were able to provide 
all non-vanishing components of the induced fermionic current densities in closed forms. The 
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corresponding results were analyzed in the massless limit and large mass limit as well. In the 
latter a strong exponential decay takes place. 

The induced fermionic current densities calculated in this paper satisfy the conservation 
condition, ^ a{j A {x)) Ren. = 0. Another important point which should be emphasized, is the 
fact that they depend crucially on the fractional part of the ratio of the magnetic flux by the 
quantum one. This is a direct consequence of the well known Aharanov-Bohm effect. 

An interesting result which deserves to be mentioned is about the calculation of charge 
density in three dimensions. As we have seen, for massless fields there is no induced charge 
density; however, considering massive fields, a non-vanishing density is obtained and an induced 
electric charge appears. This charge is a fractional part of e. 

In odd dimensional spacetime, the Clifford algebra has two inequivalent irreducible repre- 
sentations which cannot be mapped to each other by similarity transformations. They may be 
represented by El 

r = (V, 7 2 , 7 d+1 ) and j» = ( 7 \ 7 2 , - 7 d+l ) , (133) 

where j d+1 = ^/y 12 — or ^ d+1 = ±ij 12 --- d , depending on the dimensions of the spacetime, being 
7 12 -- d the product of 7*. 

In this present paper we have used 7 1+d = 7°. If we had assumed the second representation 
for the Dirac matrices, the expression found for the bispinor remains the same; however, there 
will appear a modification in the fermion propagator given in (|15p due to the presence of 7 
matrix. This new expression for Sf(x,x'), would provide a change in the sign of the charge 
density and induced electric charge for massive fermionic field case analysed in three dimensions. 

Finally we want to say that the renormalized procedure adopted to calculate the VEVs of 
fermionic current densities, was the point-splitting one. In this way we have investigated the 
behavior of the fermionic propagators in the coincidence limit and extract all the divergences in 
a manifest form by using the Haddamard functions P^l 
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A Special analisys of induced fermionic current in three dimen- 
sional cosmic string spacetime 

In this appendix we present explicitly the steps developed to prove that the radial current in 
three dimensions for massless fermionic fields, given by Eq. (155p . in fact is zero; moreover we 
show that the azimuthal current density, given for general situation by an integral representation 
in Eq. ()65|) . coincides with the (|102j) . when we substitute in the former the planar angle deficit 
and magnetic flux by specific values. 



A.l Radial current 

In order to prove that the radial fermionic current for massless fields in three dimensions is zero, 
we write down the result of the indefinite integral which appears in (|55H . It reads, 

dv 14 4cosh(t>7/a) 

-G{v, a, 7) = . , . . r) ' { ■ (134) 



sinh(?;/2) cosh(ty'a) — 1 ' ' sinh(u/2) a sinh(f/2a) 



12 See paper by Jeon-Hyuck Park in the site http://conf.kias.re.kr/~brane/2003/gamma.pdf 
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Because the cosmic string spacetime is locally flat the Hadamard function coincides with the usual Green 
function. 
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We can see that the behavior of the integral goes to zero when v — )■ and also goes to zero when 
v — > oo for 7 is smaller than ^, that is compatible with the condition assumed for the parameter 
in (|26|) . as we have already discussed in subsection 13. 1.1 1 



A. 2 Azimuthal current 

The general expression for the azimuthal current density associated with massless fermionic field 
in three dimensions can only be provided by an integral representation. In this subsection we 
shall present that if we consider the special values for the planar angle deficit and magnetic flux 
exhibited in subsection 12.21 the integral can be exactly solved and the result coincides with that 
one found in subsection 13.2.11 

First let us write ([65]) adopting a = 1/q and 7 = (q — l)/2q. The obtained expression reads: 

en 3 r°° dr 

WW) = / tttt /(*) , (135) 

QTT^r^ J_ 00 sinn(x) 

with 

f(x) = cosh(x) sinh(x) (l — coth 2 (gx)) H (coth(gx) — tanh(x)) . (136) 

In this new expression we have introduced the variable v = 2x and extend the interval of 
integration to cover all the real axis x. 

We can write f(x) as the sum of two other functions, 



with 



f(x) = f 1 (x) + f 2 (x) , (137) 

fi(x) = cosh(x) sinh(x) (l — coth 2 (gx)) H ^ — 

„ . . cosh 2 (x) . , , . , , coth(x) /^„„\ 

f 2 ( x ) = — (coth(gx) - tanh(x)) ^ . (138) 



The integral in (f!35]) . 



' Lr f(x) , (139) 



sinh(x) 



can be evaluated by using calculus of residues. 

Let us investigate the contour integral in the complex z— plane of the analytical continuation 
of the function f(x): 

r/: m . (140) 



sinh(z 



The contour to be considered is exhibit in figure 2: 

The poles of f(z)/ sinh(z) on the complex z— plane are in the imaginary axis, and are given 
by Zk = kiri/q, for k = 1, 2, q — 1. 

By the residue theorem we have 

-^-f(z) = 2^Y j a- 1 ^. (141) 
sinh(z) ^ 
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Im(z) z-plano 




Figure 2: The contour used to calculate ()140[) . 

The above integral can be evaluated by using the functions fi(z) and $2{z) separately. Let us 
start with the integral for fi(z). Because the integrals along the two vertical lines go to zero as 
R goes to infinity and for the upper horizontal line fi(z)/ sinh(z) — > — f\{z)/ sinh(z) we get: 

OO J 

^Jk) MZ)= "^ a -^- <142) 

The function fi(z) / sinh(z) presents poles of order 2. By using standard procedure to calculate 
the residue [IB], we get: 

(fc) = _ * , x / 143 v 

q z 

As to the integral involving f2(z), we also can write 

rD0 <iz 



/ 2 (z) = ^^a_!( fc ) . (144) 



5 sinh(z 

The poles of the integrand are simple ones and for this case we found 

(fe) _ i cos 2 '(kn/q) 
_1 g 2 sin(A;7r/g) 



(145) 



So taking into account the both results above for the residues of the integrals, we can see 
that 

dz 



j sinh(z) q 2 j^— J sin(kir/q) 

Substituting this result into (1135R we re obtain (I102p 
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